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■ Abstract. Gamma-ray bursts (GRBs) are a mixed class of sources con- 
(-| I sisting of, at least, the long duration and short-hard subclasses, the X- 
0^- ray flashes, and the low-luminosity GRBs. In all cases, the release of 
Q I enormous amounts of energy on a short timescale makes an energetic, 

■ relativistic or mildly relativistic fireball that expands until it reaches a 
c/3 I coasting Lorentz factor determined by the amount of baryons mixed into 

, the fireball. Radiation is produced when the blast wave interacts with 

I the surrounding medium at an external shock, or when shell collisions 

■ dissipate kinetic energy at internal shocks. This series of notes is orga- 
^ i nized as follows: (1) The observational situation of GRBs is summarized; 

(2) Progenitor models of GRBs are described; (3) An overview of the the 

, blast-wave physics used to model leptonic emissions is given; (4) GRB 

' physics is applied to hadronic acceleration and ultra-high energy cos- 

Q>^ , mic ray production; (5) Prospects for GRB physics and 7-ray astronomy 

' with the Fermi Gamma-ray Space Telescope (FGST, formerly GLAST), 

00 , and space-based and ground-based observatories are considered. Also 

included are exercises and problems. 
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Vh ' 0. Introduction 



GRBs are brief flashes of radiation at hard X-ray and soft 7-ray energies 
that display a wide variety of time histories. GRBs were first detected with the 
Vela series of spacecraft at soft 7-ray energies with wide field-of-view instruments 
used to monitor terrestrial nuclear explosions. The Burst and Transient Source 
Experiment (BATSE) on the Compton Camma Ray Observatory showed that 
GRBs are not a Galactic disk population. The discovery of X-ray afterglows 
with Beppo-SAX allowed for counterpart identification and redshift measure- 
ment, revealing the distance scale to long duration GRBs. HETE-II and Swift 
observations have revealed the counterparts to the short-hard class of GRBs, and 
Swift observations give us a crucial look at the early afterglow X-ray behavior 
of GRBs. 

This series of five lectures is organized as follows: 
1. Observations of CRBs. Introduction to GRB observations; 
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2. GRB Progenitor Models. A brief review of GRB progenitor models; 

3. Leptonic processes in GRBs-prompt and afterglow emissions. A descrip- 
tion of the blast wave physics used to model leptonic emissions, including 
nonthermal synchrotron and synchrotron self-Compton radiations; 

4. Hadronic processes and cosmic rays from GRBs. Application of blast wave 
physics to hadronic emissions and the acceleration of cosmic rays; 

5. GRB and ^-ray studies with the Fermi Gamma-ray Space Telescope, for- 
merly the Gamma-ray Large Area Space Telescope, GLAST. The future of 
GRB studies and prospects for 7-ray astronomy in light of Fermi, ground- 
based 7-ray telescopes, and complementary multiwavelength and multi- 
messenger observatories. 

The term "gamma-ray bursts" is now understood to comprise several classes 
of these sources, including 

1. Long-duration GRBs, the source class commonly meant by classical GRBs, 
which are associated with high-mass stars and star-forming galaxies; 

2. The short hard class of GRBs, which are widely thought to originate from 
the coalescence of compact objects; 

3. The X-ray flashes, which are distinguished from the long-duration GRBs 
by having the peaks of their energy output at X-ray rather than soft 7-ray 
energies; and 

4. Low luminosity GRBs, which produce unusually low energy releases com- 
pared to the long-duration GRBs. 

In addition to these subclasses of GRBs are the soft gamma repeaters (SGRs), 
which technically are not GRBs, but rather related to phenomena on highly 
magnetized neutron stars. Even without black holes they may have, like radio- 
loud subclasses of black-hole sources, relativistic outflows. 

The physics developed here can be applied to all these classes of GRBs, 
and other important cosmic phenomena (e.g., blazars, microquasars) , because 
all involve the release of a large quantity of energy during a catastrophic event, 
thought to be driven by matter accretion onto a black hole or, in the case of 
SGRs, a crustal anomaly on a highly magnetized neutron star. 

1. Gamma-ray Bursts: Overview of the Observations 

A GRB may flare up from any direction in space. The classical, long-duration 
GRB (LGRB) releases most of its energy in hard X-ray and soft 7-ray (X/7) 
energies during (to the observer) a fraction of a second to tens of seconds. There 
is no compelling evidence that LGRBs are recurrent events. Therefore, a wide 
field-of-view instrument is necessary for serendipitous detection. 
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1.1. Discovery of Gamma Ray Bursts 

GRBs were discovered in data returned between 1967 and 1973 by the Vela 
series of satellites used to monitor compliance with the nuclear test ban treaty. 
The Vela spacecraft carried non-imaging Csl detectors that were sensitive in 
the « 200 keV - 1 MeV range (Klebesadel, Strong, & Olson 1973). Above the 
large background, coincident events were identified in the light curves. Timing 
studies and triangulation were used to give an approximate direction to the 
GRBs, revealing a cosmic/non-terrestrial origin. [Exercise: Perform a simple 
timing analysis from synthetic satellite data to show how to reconstruct arrival 
direction information. Give uncertainty analysis.] 

The basic observational data in GRB studies are the spectral photon fluxes 
0(e; t) measured at time t and at photon energy hv = rrieCpe. Prom this quantity, 
one can derive, after subtracting backgound flux, the uF^ flux (cgs units of ergs 
cm~^ s~^) /e(i) = rUeC^ J^^ de e (p{e;t) « meC^e^(f){ed;t), where is the typical 
photon energy at which the detector is most sensitive. The fluence between times 
ti and t2 and within the energy range ei and 62 (cgs units of ergs cm~^) is given 
by J^{ti,t2-,ei,€2) = nieC^ J^^ dt J^^ de e (p{e;t). Flux and fluence distributions 
can be constructed from observations of many GRBs. 

1.2. BATSE Observations: GRBs are Cosmological 

The Burst and Transient Source Experiment BATSE on CGRO consisted of an 
array of large area detectors (LADs) most sensitive in the 50-300 keV band, 
in addition to smaller spectroscopy detectors. The BATSE has given the most 
extensive data base of GRB observations during the prompt phase. It searched 
for GRBs by examining strings of data for > 5.5a enhancements above back- 
ground on the 64 ms, 256 ms, and 1024 ms time scales, and triggers on GRBs 
as faint as ~ 0.5 ph cm~^ s~^, corresponding to energy flux sensitivities ^ 
ergs cm~^ s~^. At hard X-ray and soft 7-ray energies, the peak flux may reach 
hundreds of photons cm~^ in rare cases. Empirical morphological studies 
of LGRBs give various phenomenological relations, including hardness-intensity 
correlation, generic hard-to-soft evolution of epk{t), and variability-distance cor- 



Expressing the sensitivity of a high-energy radiation detector in terms of a 
threshold energy flux ^f^j. (same units as /g) imposes the condition that $ > 
^thr- For unbeamed sources with luminosity and distance c?, $ = L^/And'^, 
and the maximum source distance for a give source flux $ is 



(The luminosity distance which includes cosmological effects is defined by d^ = 
y'L*/47r$, and d = at low redshifts z <^ 1. [Exercise. Relate energy to 
fluence, including redshift.]) Hence the well-known —3/2 result for sources uni- 
formly distributed with density uq in Euclidean space, namely 



relation. 




N{> $) = N{< d) = 47rno / 



.d(<i.) 
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The {V/Vmax) statistic 

expresses the deviation from 0.5 expected for a uniform Euclidean distribu- 
tion of sources (Schmidt 1968). Here V stands for volume, and Vmax is the 
maximum volume from which a source with flux $ could be detected. Val- 
ues of (V/Vmax) > 0.5 represent positive evolution of sources, that is, either 
more sources and/or brighter sources at large distances or earlier times. Val- 
ues of {V/Vmax) < 0.5 represent negative source evolution, i.e., fewer or dimmer 
sources in the past. Detailed treatments of the statistical properties of black hole 
sources must consider cosmological effects and evolution of source properties. 

The integral size distribution of BATSE GRBs in terms of peak flux (pp is 
very flat below ~ 3 ph cm~^ s~^, and becomes steeper than the —3/2 behav- 
ior expected from a Euclidean distribution of sources at (f>p ^ 10 ph cm~^ s~^. 
The directions to the BATSE GRBs arc isotropically distributed in the sky and 
display no clustering toward the Galactic plane. When coupled with the flat- 
tening of the peak flux distribution, the implication is that we are at the center 
of an isotropic though bounded distribution of GRB sources. A cosmological 
distribution of sources is most compatible with these observations. 

The duration of a GRB is defined by the time during which the middle 50% 
(tso), 90% (tgo), or i% ((ij)) of the counts above background are measured. A 
bimodal duration distribution is measured, irrespective of whether the ^50 or ^99 
durations arc considered (Kouvcliotou et al. 1993) About two-thirds of BATSE 
GRBs are long-duration GRBs with tgo ^ 2 s, with the remainder comprising the 
short duration GRBs. The short duration GRBs tend to have harder spectra, so 
that they are referred to as the short, hard class of GRBs. They are also much 
weaker in average fluence in the BATSE range. 

GRBs typically show a very hard spectrum in the hard X-ray to soft 7-ray 
regime, with a photon index breaking from ~ — 1 at photon energies Eph ^ 50 
keV to a —2 to —3 spectrum at Eph ^ several hundred keV. In BATSE studies, 
the distribution of the peak photon energies Epk of the time-averaged i^Fi, spectra 
of BATSE GRBs are typically found in the 100 keV - several MeV range. The 
time-averaged or, for very bright GRBs, time-sliced GRB spectrum is usually 
well-described by the "Band function" Nsie) (Band et al. 1993), a power-law 
times an exponential that smoothly connects to a steeper power-law, given by 

N^{e) = k, e°exp[-6(a - (3)/eJ H{e-el,^,eJ 

+ k, e^-^ expifi a) H{e; etr, ^Lj • (3) 
Here a and j3 are the low and high energy photon number indices, and £^br = 
me(?e^^ is the "break energy." The Heaviside function H{x;y,z) vanishes ev- 
erywhere except at y < X < z, where it equals unity. The term kg is the 
constant normalizing the number fluence to the > 20 keV BATSE energy flu- 
ence ^g{> 20 keV) of a particular GRB [Exercise: Derive the form of the Band 
function and normalizing constant, and convolve with nontrivial model detector 
response.] Typical values of Band alphas a « — 1 and Band betas /? « —2.2 

2.5. Deviations of these values give valuable information about radiation 

processes and existence of separate radiative components. 
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1.3. Beppo-SAX and the Afterglow Revolution 

Beppo-SAX GRB observations reveal that essentially all long-duration GRBs 
have fading X-ray afterglows. Beppo-SAX, launched April 30, 1996, carried 
three instruments. The Gamma Ray Burst Monitor was sensitive in the range 
60 - 600 keV to GRBs brighter than f» lO''^ ergs cm^^ g-i ^he Wide Field 
Camera was sensitive in the range 2-30 keV down to ~ 10"^'^ ergs cm~^ 
and provided ;$ 10' error boxes. The spacecraft was then slewed, requiring 6 
- 8 hours, but was fast enough for the Narrow Field Instruments, sensitive to 
> 0.1 keV GRB emissions as faint as ~ 10~^^ ergs cm~^ s~^, to give error 
boxes ;$ 0.5'. The first X-ray afterglow was obtained from GRB 970228 (Costa 
1997), which revealed an X-ray source that decayed according to a power-law, 
(j)x{t) oc t^, with X ~ —1.33. Typically, x ~ —1.1 to —1.5 in Beppo-SAX ~ 2 - 
10 keV X-ray studies. 

The small X-ray error boxes allowed deep optical and radio follow-up stud- 
ies. GRB 970228 was the first GRB from which an optical counterpart was ob- 
served (van Paradijs et al. 1997), and GRB 970508 was the first GRB for which 
a redshift was measured. Redshifts are provided by detection of optical emission 
lines from the host galaxy and absorption lines in the fading optical afterglow 
due to the presence of intervening gas. Host galaxies of long duration GRBs are 
bluish star-forming galaxies, primarily consisting of dwarf irregular and spiral 
galaxies. No optical counterparts were detected from approximately one-half of 
Beppo-Sax GRBs with well-localized X-ray afterglows, and are termed "dark" 
bursts. These sources may be undetected in the optical band because of dusty 
media, or intrinsically faint afterglows. These results give compelling evidence 
that LGRBs are associated with star-forming galaxies and the deaths of mas- 
sive stars, especially given the detection of supernova emissions a few weeks 
after the GRB in a few, nearby faint GRBs which, however, may not be fully 
representative of the LGRBs (but rather the low luminosity GRBs, LLGRBs. 

Apparent isotropic energy releases of LGRBs arc enormous, exceeding 10^^ 
ergs, and the redshift distribution of Beppo-SAX, BATSE, HETE-II and IN- 
TEGRAL (pre-Swift) GRBs is peaked near {z) ~ 1. Achromatic breaks in the 
optical curves of GRB afterglows gives evidence for a beamed/jetted geometries, 
reducing the apparent energy release to a beaming-corrected energy release by 
a beaming factor J^bm- Approximately 40% of GRBs have radio counterparts, 
and the transition from a scintillating to smooth behavior in the radio afterglow 
of GRB 980425 provides evidence for an expanding source. For BATSE/Beppo- 
SAX type GRBs, the lion's share ~ 65% of the energy is released in the form 
of > 25 keV X-rays and soft 7 rays, ~ 7% is softer X-rays, ~ 0.1% in optical, 
during the prompt phase t ^ 2{ti). 

A class of X-ray rich GRBs, with durations on the order of seconds to 
minutes and X-ray fluxes in the range 10~^ - 10~^ ergs cm~^ s~^ in the 2-25 
keV band, was detected with many X-ray satellites, including Ariel V, HEAO-1, 
ROSAT, and Ginga, but conclusively established with Beppo-SAX (Heise et al. 
2001). These X-ray flashes (XRFs) are "7-ray challenged," as indicated by Band- 
function model fits to the prompt emission spectrum. Several phenomenological 
correlations of XRFs and LGRBs have been reported. One goal is to establish a 
pseudo-redshift indicator, another a "pulse" paradigm, and correlations between 
the duration of quiescent and subsequent pulse periods in separated pulses. A 
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quantitative relation between integrated fluence and Epk in well-defined pulses 
is reported. 

The Amati relation (Amati et al. 2002) correlates the i^Fi, peak photon 
energy Epk with apparent isotropic energy release Siso = 10^^£^54 ergs according 
to 

Epk oc . (4) 

The Ghirlanda relation (Ghirlanda et al. 2004) correlates Epk with the collimation- 
corrected absolute X/7 energy release Sabs = 10^-^ £^51 ergs according to 

Epk oc ■ (5) 

The Amati and Ghirlanda relations are challenging to explain and potentially 
useful for cosmolgical studies. 

1.4. Swift and Various Classes of Bursting Sources 

The Swift Observatory is a NASA-ASI-supported MidEx launched November 
20, 2004. Its main scientific payload consists of 3 instruments: the Burst Alert 
Telescope (BAT), triggering between 15 and 150 keV, and a very sensitive X- 
ray Telescope (XRT) operating between ^ 0.3 - few keV, and the Ultra-violet 
optical telescope (UVOT). The spacecraft slews autonomously in 20 - 75 s in 
response to triggers from the BAT. 

Long duration GRBs (LGRB). The redshift distribution of the GRBs detected 
with Swift, with average Swift GRB redshift {z) 2, differing markedly from 
(z) ~ 1 for prc-Swift GRBs. This can be explained by the difi'crcnt triggering 
energy range of BAT vs. BATSE, and Epk-Hux correlations. Knowledge of the 
early X-ray afterglow phase, an important goal realized by the Swift mission, 
provided surprising unpredicted behavior, most notably rapid X-ray declines 
and X-ray flares to late times (Zhang et al. 2006; Nousek et al. 2006). The 
physical meaning of the rapid declines in the X-ray flux between ~ 100 - 10^ s 
is in dispute, including muti-component jet models, refreshed jets, and hadronic 
signatures. 

Short Hard GRBs (SGRB). As expected for an old population of host galaxy 
progenitors, as would be the case for coalescing neutron stars and black holes, 
SGRBs are expected to be associated with ellipticals as well as spirals. No ev- 
idence for supernova emissions has been found in SGRB afterglows, meaning 
that the progenitors are not associated with a young, massive stellar popula- 
tion. That this expectation received spectacular confirmation with Swift, and 
also HETE-II, only opens up the surprises, including delayed X-ray afterglow 
emissions, rapid X-ray declines, and X-ray flares. 

The heterogeneous Swift/HETE sample has a redshift distribution broadly 
distributed around {z) « 0.4, and differs in important ways from the LGRBs 
in terms of host galaxies, ofi'sets from host galaxies, apparent energy releases, 
and lag-luminosity relation. Absolute energy determination for SGRBs is com- 
promised by the difficulty of finding beaming breaks in light curves of SGRBs. 
(Exercise: Beppo-SAX was not sensitive to SGRBs. Why?) 
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Low Luminosity GRBs (LLGRBs). The LLGRB class is the long tail to the 
LGRBs to low apparent energy release (reviewed by Zhang 2007). LLGRBs 
could very well be a separate and distinct population from LGRBs, possibly 
associated with magnctar activity, as indicated by GRB 060218, or extended 
black-hole fueling activity. A possible new class of GRBs was found in relation 
to nearby GRB 060614 with no supernova emissions, though it could possibly be 
a nearby SGRB. Identification of LLGRBs as a separate class, known since GRB 
980425/SN 1998bw, makes us rethink conclusions about the relation of LGRBs 
to supernovae (SNe), which mostly derived from nearby LLGRBs. {Exercise: 
Report on GRB 98025.) 

Soft Gamma Repeaters (SGRs). The giant fiare of December 27, 2004 was the 
most intense cosmic transient observed historically. It was detected by over 20 
spacecraft from the Earth to Saturn, and apparently released ~ 1000 times more 
energy than all the Milky Way's stars, ^ 1 erg cm~^ in hard X-rays and 7 rays 
at Earth (Hurley et al. 2005). {Exercise. Estimate the bolometric radiation flux 
from the stars in the Milky Way. Estimate the total stellar energy flux of the 
universe.) The event, releasing ~ 10^^ ergs, lowered the level of the Earth's 
ionosphere. 

SGR 041227 began with an ~ 0.2 s long, hard spectrum spikes with E ~ 
-[^q46 _ -[^q47 gj.g_ 'pj^g spike is followed by a pulsating tail with ~ I/IOOO*'' of the 
energy. Viewed from a large distance, only the initial spike would be visible, 
and would resemble a short GRB. It could be detected out to 100 Mpc GRB 
050906 at z = 0.03 could be a magnetar flare. Giant flares like this must occur 
in other galaxies, and could comprise ~ 10% of the the SGRB population. 

1.5. Summary 

1. Pioneering phase (1967 - 1991): era of confusion 

2. BATSE/CGRO era (1991 -2000): GRBs are cosmological 

3. Beppo-SAX afterglow era (1996 - 2006): distance and energy scale for 
classical LGRBs (HETE-II/INTEGRAL) 

4. XRFs, first seen with Ginga, clarified with Beppo-SAX 

5. Low luminosity GRBs class (first example in 1998: GRB 980425/SN 1998bw) 

6. Swift era (2004 - ): early afterglows and the distance and energy scale for 
short-hard SGRBs (HETE-II) 

7. Fermi GBM (GLAST Burst Monitor) /Fermi LAT (Large Area Telescope) 
era (2008 - ) 

2. GRB Progenitor Models 

Summarizing the observations of LGRBs, we find that a "typical" long-duration 
GRB lasts s» 20 s in hard X-ray/soft 7 ray emission from keV to MeV energies. 
It takes place in star-forming (spiral or dwarf irregular) galaxies, but not in 
ellipticals. It takes place in a galaxy at (z) 1 - 2, and releases ~ 10^^ - 10^^ 
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ergs of apparent isotropic energy in bursts of radiation with apparent isotropic 
luminosities of « 10^° - 10^^ ergs/s. It is followed by long-lived X-ray, optical, 
and radio afterglow emission. Variability times arc as short as ms (though more 
typically 1 s) . How to explain this phenomenology? The consensus of most GRB 
scientists is that a LGRB is the consequence of the collapse of a massive star 
Ri SOMq to a black hole fed by an accretion torus. 

In this second Lecture, we present a description of 

1. GRB Source Models: Core Collapse vs. Coalescence 

2. Energy source: nuclear, gravitational, rotation 

3. Jet formation: neutrinos vs. magnetic field energy 

4. Pathways to progenitor formation of collapsars 

5. Pathways to progenitor formation of merging compact objects 

2.1. Beaming and the Fireball Model 

A strongly collimated, jetted emission helps answer the question, "How is it 
possible for a source to produce 10^^ ergs of energy?" The answer is, you produce 
only a fraction of this energy, in beams with dE/dO, E/A-k in a small solid 
angle element AO. Beaming into a cone of opening half- angle 0j « 5 - 10° 
decreases the gamma-ray energies by ~ 2 orders of magnitude, to ~ 10^^ erg. 
[Exercise: Work out the beaming factor for a two-sided top-hat jet, and for 
a non-trivial jet profile.] Optical beaming breaks are derived by equating the 
characteristic beaming angle of emission, s» l/F, with the jet opening angle, 9j. 
A softening of the power law takes place when Oj ^ l/F as deceleration of the 
blast wave reduces the T factor. Beaming also increases the total burst rate 
by the same factor, but this does not contradict anything we know about star 
formation and evolution for J~\)rn ^ 

10^ - 10^. 

But how to obtain large energy releases, large luminosities, and short vari- 
ability time scales?, and 

How does the kinetic energy of the ejecta get converted to electromagnetic 
radiation? 

These questions are answered by the fireball/blast wave model considered 
in the next lecture. In this lecture, we instead consider the question of the 
progenitor to the 7-ray burst. 

Compared to SNe, which occur ~ once per sec throughout the universe, 
LGRBs take place, depending on beaming factor, ~ J'}mi P^r day « 102(^6^/100) 
per day throughout the universe, representing as few as 0.1% of the SN popu- 
lation. GRB-hosted SNe, for example, SN 2002ap, tend to by hyper-energetic 
compared to "normal" SNe. Thus Paczynski (1998) coined the name "hyper- 
novae," which has proven to be a useful concept. 

2.2. Classes of 7 Ray Transients 

Short GRBs. What are the mechanisms for the SGRBs? Most GRB scientists 
support a picture involving the merger of a compact binary system, such as 
two neutron stars, or a neutron star and a black hole (Eichler et al. 1989). 



La Plata Lectures on GRBs and Fermi/GLAST 



9 



Lack of optical and radio afterglows is explained by tenuous ISM, if the merger 
takes place outside the host galaxy. The class of GRBs involving merger events, 
namely the SGRBs, arc different in significant ways, clearly revealing the two 
as having distinct origins. But there is also another distinct class. 

X-Ray Flashes. [Exercise: Report on the properties of the XRFs.] 
If XRFs are another manifestation of long GRBs, then are they 

1. GRBs at high redshift? 

2. GRBs observed away from the jet axis? 

3. Explosions with less relativistic ejecta? 

We have fairly complete data on one XRF (XRF020903, z = 0.251); in this 
case, the answer is compatible with 3 (Soderberg et al. 2004). The Amati relation 
gives evidence that XRFs and GRBs are part of the same family, consistent with 
a mass-loaded fireball model. 

Soft Gamma Repeaters. What are the SGRs? Most GRB scientists think that 
the 3 or 4 SGRs in the Milky Way and the SGR in the LMC originate from 
a transient release of magnetic field energy in highly magnetized neutron stars 
{B > Bcr = mlc^/eh). The energy could be released from global crustal fracture 
through B field annihilation, or when the magnetosphere fills with hot e~e"'" 
plasma. 

Long duration GRBs. LGRB progenitors generally involve black- hole/accrction- 
disk models. To answer the energy generation and release question, a number 
of problems must be solved, including the problems of the 

• Black Hole/ Accretion disk (BHAD) + Jets; 

• Structure of a Relativistic Disk; 

• Neutrino Driven Explosions from a BHAD system; and 

• Magnetic Field Driven Explosions from a BHAD system. 

To firmly resolve the astrophysical origins of GRBs, it is necessary to also under- 
stand the connections of the various classes of GRBs to their stellar evolutionary 
past. 

2.3. Energy Sources and Energy Transport 

Three main energy sources in astrophysics are 

1. Nuclear Energy 

2. Gravitational Potential Energy 

3. Black Hole Rotational Energy 
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The proton-proton chain efficiency is 0.0067, too small to reahstically power 
most GRB radiations. Gravitational energy will be released during core collapse 
or coalescence event, and such events could also trigger a rapid discharge of 
rotational energy from a rapidly spinning black hole. Which of the other two 
are most important for GRB energy release, or are they both? 

Energy in the form of photons is transported very differently through stellar 
interiors than the energy of neutrinos, the former dominated by absorption and 
reemission, and the latter often following a rectilinear trajectory out from the 
stellar core. Simulations of core collapse give us ways to follow different modes 
of energy release, also including gravitational radiation. 

2.4. Principal GRB Models 

The majority GRB establishment has settled on a standard scenario for the pro- 
duction of LGRBs and SGRBs, both involving black-hole accretion disk (BHAD) 
models. The energy release in LGRBs is believed to result from the accretion of 
a massive torus onto a black hole following the collapse of the core of a massive 
star to a black hole. This is the collapsar scenario, pioneered by Stan Woosley 
(1993). Compact object mergers have a different family history, but also end up 
with the release of energy by jetted outflows from BHAD-type systems. 

A variety of scenarios for catrastophic events involving compact stars can 
be imagined, for example, 

1. Neutron Star - Neutron star mergers, e.g., the Hulse- Taylor binary pulsar 
system 

2. Black Hole - Neutron star mergers 

3. Black Hole - White dwarf mergers 

4. Collapse of the core of a rotating massive star (binary or single star) 

5. Neutron star - black hole merger with the helium core of an evolved pri- 
mary, a He-star merger 

The collapse of the core of a rotating massive star down to a black hole is 
the original scenario behind the collapsar model. MacFadyen & Woosley (1999) 
modeled the toy case of such a collapse, placing an artificial angular momentum 
profile onto a massive star whose core was assumed to immediately collapse to 
a black hole. They identified many of the key ingredients for a succesful GRB 
outburst: high accretion rate and specific angular momentum. This simulation 
provided the first numerical confirmation of the collapsar model and showed that 
with artificial, but physically reasonable, energy drivers that the collapsar could 
explain the durations of LGRBs. It is worth noting that since this calculation, 
much more physically relevant models have been run including either magnetic 
fields (Proga et al. 2003) or modeling the collapse in 3-dimensions of a realistic 
star (Rockefeller et al. 2006). 

One way to drive the energy is through neutrino annihilation above the 
black hole accretion disk. The usual assumed geometry of a torus accreting 
onto a black hole is given by a wedge, with prolific neutrino fluxes that forms a 
relativistic e''"e~ pair wind from neutrino-antineutrino interactions. In extreme 
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cases for a large Shakura-Sunyaev accretion disk a(« 0.1) parameter and a 
maximally rotating black hole, efficiencies for i/v' — >e"'"e~ could reach a few 
percent. This channel for energy release is thought to be too inefficient to power 
LGRBs, though could be important for subclasses of SGRBs. Large a, rapidly 
rotating black holes, can release sufficient energy to power a GRB. geometric 
beaming can be produced from energy deposition via neutrino annihilation. See 
Popham et al. (1999) for a review. 

As far as neutrino-driven GRB engines go, the critical densities for most 
likely accretion disks are ~ 10^ - 10^ g/cm^. For collapsars, this corresponds 
to black hole masses of « 10 - 25Mq, with delays between collapse and jet 
formation of 30 - 300 s. The neutrino-driven collapsar model probably does not 
work (Fryer & Meszaros 2003). The alternatives arc magnetic fields and fast 
rotation (Narayan et al. 1992). Again, see Popham et al. (1999) for a review. 

The magnetic field in the accretion-disk/black-hole/jet system is probably 
dominated in currents flowing in the accretion disk. Penrose energy extraction 
through the Blandford-Znajek (BZ) process can be as powerful as the accre- 
tion power. The literature on magnetically driven jets is extensive. The BZ 
power is insufficient except for rapidly rotating black holes. Magnetically driven 
jets probably produce much more energy than neutrino annihilation. Common 
physics should apply to GRBs and AGNs with collimated relativistic outflows, 
namely radio galaxies and 7-ray blazars. 

Numerical simulations of collapsar jets through stellar envelopes of massive 
progenitor stars show the need for the envelope to have a width ^ 10^^ cm. 
Emerging jets, possibly protected by a cocoon formed at the interface of the hot 
light relativistic fluid, can maintain relativistic speeds. In the picture of Woosley, 
LLGRBs such as GRB 980425 is an off axis GRB, and XRFs could arise from an- 
gular effects in a mildly relativistic fireball, in contrast to a "dirty," mass-loaded 
fireball scenario by Dermer, Bottcher and Chiang (1999). (Exercise: Tabulate 
similarities and differences, and contrast predictions of off-axis vs. mass- loaded 
XRF scenarios.) Constraints on the formation of collapsar GRBs include: 

• The star must collapse to form black hole. 

• Star must lose its hydrogen envelope so that it remains compact. 

• Jet must travel through star roughly on the GRB duration timescale. 

• Star must be rapidly rotating so disk forms around black hole. 

An important mechanism in stellar evolution leading to GRBs is Roche- 
lobe overflow. When a star expands (in the giant or supergiant branch) in a 
binary, the outer layers of the star may feel greater gravitational attraction to 
the companion star, causing this material to "overflow" onto that companion. 
A second important mechanism is common envelope (CE) evolution. If Roche- 
lobe overflow proceeds faster than the companion star can accrete material, the 
expanding star envelops the companion, leading to two stellar cores within a 
single common envelope. The masses of primary (most massive) and secondary 
(least massive) stars in a binary are Mp and Mg. 

Some pathways to GRB formation include: (i) Massive Star which, if not 
affected by binary mass transfer, would undergo core-collapse {Msn ~ 8 - 10 
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Mq) to a black hole; (ii) He core - He core mergers during CE phase of binary 
systems with two massive stars (He core masses will also have transition masses 
for neutron star and black hole formation). The Black Hole Mass {Mbh ~ 
3M0) is the transition mass for black hole formation, and depends on angular 
momentum of core. See Fryer, Woosley, and Hartniann (1999). 

In the isolated collapsar scenario, a massive single star with high metallicity 
loses its hydrogen envelope via winds. If it retains enough mass and rotation 
to form a BHAD system, a GRB is produced. Important related issues are 
the metallicity of the host environment, and whether the progenitor has the 
required rotation. Modeling the full collapse (first to a neutron star and then 
to a black hole) is required to understand some of these trends. Stellar collapse 
with nuclear equations of state are now being done in 3-dimcnsions (e.g. Fryer 
and Warren 2002) including rotation (Fryer and Warren 2004) and finally to 
the late times that approach black hole formation (Fryer and Young 2007). 
From these we can trace fluid flow and isopressure surfaces of outward moving 
bubbles. {Exercise: Use the concepts of mixing length in convection to write 
the one-dimensional stellar structure equations for a collapsing object. Discuss 
the equations' deficiencies.) 

The combined simulation and theory understanding of stellar collapse ar- 
gues that different progenitor scenarios lead to remnant black holes with a mass 
distribution possibly reflected in systems with inferred black hole masses (Fryer 
& Kalogera 2001). Because large a black holes allow greater energy extraction 
per unit mass, GRB-hosted SNe are bound to be more asymmetric than normal 
SNe to have a spun-up black hole remnant, though all SNe should exhibit asym- 
metries. Note that those objects that do form GRBs eject the star, producing a 
lower mass black hole remnant than those SNe and failed SNe that form black 
holes with no subsequent explosion (presumably due to slow rotation speeds). 

One of the observed properties of GRBs is that the associated supernovae 
do not show any hydrogen. Strong stellar winds can remove the outer hydrogen 
envelope, but these winds also carry away angular momentum. An alternative 
is to assume that the star is in a binary and a common envelope phase ejects the 
hydrogen envelope. This evolution can even increase the angular momentum. 
The mass and rotation constraints are easier to satisfy in binary systems (Fryer 
et al. 1999). 

enough angular momentum. An alternative, proposed by Fryer et al. (1999) 
was the merger of two stars. Fryer & Heger (2005) merged two helium cores and 
found that they could indeed increase the core rotation rate while removing the 
hydrogen and some of the helium envelopes. For collapsars, this means that 
with magnetic fields (or with the current prescriptions of magnetic fields) , single 
stars and even simple binaries have trouble working (Petrovic et al. 2005). Yoon 
&L Langer (2005) have argued that single stars at low metallicity can go through 
extensive mixing, removing all of the hydrogen by burning it into helium instead 
of through a wind. 

But any binary may well not work. The issue is that we not only don't see 
any hydrogen, but we don't see helium in the supernovae associated with GRBs. 
So the models still have some work before they can fit all the data (see Fryer et 
al. 2007 for a review). 
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By comparison, compact object coalescence scenarios require a binary (or 
multiple star) system where the primaries and secondaries evolve to compact 
objects and merge within a Hubble time. 

Double neutron star formation for SGRBs involve neutron star formation 
where the primary collapses after expanding off the main sequence, or being 
triggered by Roche lobe overflow of material onto the massive core. A likely 
scenario in close binaries involves CE evolution until the compact object settles 
to the center of the companion. The merger of the compact object with a He core 
forms a GRB. Alternately, in the absence of CE evolution, a compact binary can 
form with two progenitor high-mass stars. Accretion-induced collapse during the 
CE phase or during Roche lobe overflow introduces other evolutionary pathways 
to CRBs. {Exercise. Report on a GRB progenitor pathway in more detail.) 

Optimistically, a neutrino mechanism could power a GRB, though would 
be associated with an old stellar population and large galactic offsets. Fryer et 
al. (1999) concluded that for collapsars, it is difficult to make isolated stars with 
enough angular momentum after removing envelope with winds. Binaries seem 
necessary (e.g. He-He or He-neutron star mergers). Compact object mergers 
(NS/NS or NS/BH) are rare, but sufficient for GRBs. Accretion rates are high 
( > 1 - IOMq/s), but short-lived ( < 200 ms). 

GRB studies intersect many exciting fields of research, including 

1. Jet and fireball physics; 

2. Ultra-high energy cosmic ray acceleration and neutrinos; 

3. Early universe, reionization; 

4. Gravitational radiation; and 

5. Mass extinctions, and geophysics. 
Competing models include the 

1. Supranova and two-step (SN ^ NS ^ BH) collapse processes for LGRBs, 
originally proposed by Vietri and Stella (1998,1999); see Dermer (2008); 
and 

2. AIC of neutron stars to black holes in binary systems as a model for SGRBs 
(MacFadyen et al. 2005; Dermer and Atoyan 2006). 

3. Leptonic processes in GRBs - prompt and afterglow emissions 

This lecture gives an overview of the elementary leptonic blast wave physics 
developed to explain GRB prompt and afterglow emissions. The outline of this 
lecture is 

1. Relativistic Kinematics 

2. Blast- wave and Afterglow Theory 

3. Relativistic Shock Hydrodynamics 
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4. Jetted Emission and Beaming Breaks 

5. External Shocks 

6. Colliding Shells 

Reviews of GRB blast wave physics are given by Meszaaros (2006), Piran (2005), 
and Dermer & Menon (2009). 
To recap, we have associated 

1. LGRBs Collapsars 

2. SGRBs ^ Mergers of compact objects 

3. LLGRBs Probably a type of collapsar, perhaps magnetar-powered 

4. SGRs ^ Highly magnetized neutron stars 

GRBs must be Galactic, as we now show based on discoveries made in the 
1980s (and stih valid). The typical energy flux of a GRB is $ = 10"^$_6 ergs 
cm~^, with significant factors-of-2 variations on timescales of seconds. Solar 
Maximum Mission observations of MeV photons (Matz ct al. 1985) showed that 
the optical depth for attenuation of 7 rays with energies ^ 1 MeV is t^^ « 
n^a^^R ~ rijaTR < 1 for photons above threshold ee' 2. The photon energy 
density ~ u^/Ej, ~ L^tesc/V ~ 3(P^/R^c, with escape time tesc ~ R/c. 
Hence 

c-'AtE^ ' 

or d{10 kpc)< ^At(s)/$_6. GRBs are known to be at cosmological distances. 
The flaw in this argument is the assumption that the emitting region is at rest. 
(The lesson here is not "Don't believe theorists.") 

The conventional fireball/blast-wave model involves intermittent release of 
energy in a collimated relativistic jet, with variations of wind parameters lead- 
ing to internal shocks and structure in the light curves. External shocks are 
responsible for the late Beppo-SAX afterglow, and could contribute to emission 
in the prompt and early afterglow phase. 

The central problem in GRB astrophysics is to explain large apparent 
isotropic energies ^ 10^^ ergs, durations ^ 10^ s, and short variability time 
scales At ;S 1 s. The widely accepted solution is the fireball/relativistic blast- 
wave model. The impulsive release of a huge amount of energy in a fireball with 
large entropy per baryon is described by an 

1. Expansion phase; T{x) f« x/Ao,x < FqAo; a 

2. Coasting phase; T{x) = Tq,x > FqAo; and a 

3. Deceleration phase; r(x) =? 
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3.1. Relativistic Kinematics 

Important questions include 

• How to calculate r(x)? 

• How is r(x) related to observer time t? 

• How to calculate internal photon number and energy density for a rela- 
tivistically moving source from measured energy flux? 

We consider these problems for a simple blast-wave model with several simpli- 
fying assumptions: a spherical, uncollimated explosion, a uniform surrounding 
medium, the blast wave approximated by a uniform thin shell, and particle ac- 
celeration at the forward shock only. Three frames of reference are considered: 
the explosion (GRB, stationary, or starred) frame; the comoving fluid (primed) 
frame; and the observer (unstarred) frame. 
Using the definition of the Doppler factor 

du = [r{l-(3ii)]-\ (6) 

where ji = cos 6 and 6 is the angle between the direction of the radiating fluid 
and the observer, we have from elementary considerations 

dt = {l + z)^, and e=p^. (7) 

For well-defined pulses of radiation on a measured variability time scale tyari 

the 

comoving emission-region size scale 

The variability timescale t^ar implies an engine size scale, with the comoving 
size scale a factor ~ Sy) larger, and the emission location ^ larger than 
values inferred for a stationary region. Rapid variability can be obtained from 
an external shock or internal shock by energizing regions within the Doppler 
cone, as defined by the Doppler factor of the shocked fluid. 

The vFi, spectrum at dimensionless photon energy e = hv/nieCp from a 
source at redshift z and luminosity distance (ii,(z) = 10^^(i28 cm is denoted by 
/g. Assuming that the blob radiates isotropically in the comoving frame. 



f (t) ^ e'i'(e'-t') ^ ^D^'^'i^''^^'^ (Q) 



where the comoving spectral luminosity L'{e';t') = V^j'{e';t'), and j'{e'\t') is 
the comoving spectral emissivity. The spectral energy density u'(e') = t'iJ'{e'). 
Because Vl = Attt'^/S, 

u', = e'u'ie') = m.cW{e') ^ |^ ^ ^fM^^ " ^''^ 
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Consequently the relation between the total internal photon energy density u' 
and the total measured energy flux <I> is 

> 34(1 + ^)^^ (11) 

Requiring that r^^(ei) < 1, so that the emission region is transparent to 
7 rays gives, using a 5-function result for the 77 attenuation cross section, the 
result 

cj'P 2,2^ 

Tjyi^i) ~ -:ri—)'^ph{—yb ■ (12) 
J ^1 ^1 

With eq. ([8]) , the requirement of optical thinness to 77 attenuation gives for flat 
target photon SEDs in a I'Fiy representation the result 

6^> 200[{l + z)d,sf'[^^^]'/'. (13) 

important for interpretation of Fermi Gamma ray Space Telescope (FGST) re- 
sults. 

3.2. Blast Wave and Afterglow Theory 

The GRB explosion forms a radiation-dominated fireball with injection explosion 
entropy per baryon % = L/Mc^ ^ 1, and L is the wind power. The energy of 
the expanding relativistic wind is transformed into photospheric emission and 
the directed kinetic energy of a hadronic shell with coasting Lorentz factor 

To = Eo/Moc^ . (14) 

Here Mq is the amount of baryonic matter mixed into the initial explosion. For 
a uniform spherically symmetric CBM, the mass of swept-up material at radius 
X is Msw = 47rmpnox^/3, where no is the proton density, assumed to be made 
of H. 

The blast wave will start to undergo significant deceleration when an amount 
of energy comparable to the initial baryon energy Eq in the blast wave is swept 
up. Looked at from the comoving frame, each proton from the CBM carries with 
it an amount of energy TorripC^ when captured by the blast wave. After capture 
and isotropization, the amount of energy carried by the blast wave from this 
swept-up proton is T^mpC^ as measured in the stationary frame. The condition 
TqMswC^ = Eq gives the deceleration radius (Rees & Meszaros 1992; Meszaros 
& Rees 1993) 

X. . ( ^ )V3 ^ 2.6 X 10^^(-^)V3 em , (15) 

where £^0 = £^52/10^^ ergs is the total explosion energy including rest mass 
energy, Faoo = ro/300, and no is the CBM proton density in units of cm"^. 

Differential time elements in the stationary (starred), comoving (primed), 
and observer (unscripted) reference frames satisfy the relations 

dx = Pcdt^ = pVcdt' = (3c ^ 



{l + z){l-(3^i) 
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where the last expression is obtained by noting that dt/{l + z) = dt'{l — jdfi), 
and 9 = arccos n is the angle between the direction of outflow and the observer. 
Hence 

dt=^-^dxip-'-^)^^^. (16) 

The last expression applies to rclativistic flows (T ^ 1) observed on-axis, as- 
suming that the average emitting region is at = /?. 

The deceleration time as measured by an observer is therefore 

3.3. Blast- Wave Equation of Motion 

The equation describing the speed of the relativistic blast wave, which changes 
as a consequence of the blast wave sweeping up material from the surrounding 
medium, is for an adiabatic blast wave, 

r[Mo + rm{x)] ^ r[Mo + ^^^(r - l)] ^ const , 

where m{x) is the swept-up mass. For r(a;) oc x~^^'^, t = f dx oc 
/ dx x^, so x{t) oc t^/^ and T{t) oc t^'^^^. If the blast wave is partially or highly 
radiative, different behaviors follow. [Exercise: Derive the power-law behavior 
for adiabatic blast waves decelerating in an external medium with radial wind 
density profiles.] 

The kinetic energy swept into the comoving fluid frame per unit proper time 
at the forward shock is given by 

J pi/ 

— \FS = A{x)nompC^{Pc)T{T-l) oc F^ for r»l. (18) 

where the area A{x) = Attx^ for an isotropic blast wave. The factor F represents 
the increase of external medium density due to length contraction, the factor 
(F — l)m,p is the kinetic energy of the swept-up particles, and the factor /3c is 
proportional to the rate at which the particle energy is swept into by the blast 
wave. This process provides internal energy available to be dissipated in the 
blast wave. The original treatment of adiabatic and radiative relativistic blast 
waves using a fluid dynamical approach was given by Blandford and McKee 
(1976). 

A fraction of the forward-shock power is assumed to be transferred to 
the electrons, so that 

(19) 

If all the swept-up electrons are accelerated, then joint normalization to power 
and number gives 

Imin = ee(^)(^)(r - 1) = 6e(^)(^)F , for F » 1 (20) 

and 2 < p < 3. The magnetic-field energy density ub = B'^/Stt is assumed to 
be a fixed fraction eg of the downstream energy density of the shocked fluid. 
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Thus B'^/Stt = AeBnompC^r^ . A break is formed in the electron spectrum at 
cooling electron Lorentz factor 7c, which is found by balancing the synchrotron 
loss time scale t'^y^^ with the adiabatic expansion time t'^^^ = x /Tc = Tt = t'gy^ = 
(4c(TT-B^7c/247rmeC^)~^, giving 

- Tfi r'^. ■ (21) 

IbeBnonipCdTi t 

For an adiabatic blast wave, T oc t~^/^, so that oc t~^l^ and 7c oc t^/^. 

The observed vFi, synchrotron spectrum from a GRB depends on the ge- 
ometry of the outflow. If Lgy^{e') = e'{dN' /de'dt') is the spectral luminosity 
in the comoving frame, then e'L'^,y„(e') = ^UBcaT^^Ne{'y), with 7 = ^e' /es 
and €b = B/Bcr = .6/(4.41 x 10^"^ G). For a spherical blast-wave geometry, 
the spectral power is amplified by two powers of the Doppler factor 5 for the 
transformed energy and time. The vF^ synchrotron spectrum is therefore 



For a power-law injection spectrum, the cooling comoving nonthermal elec- 
tron spectrum can be approximated by 

In the slow cooling regime, s = p, 7o = 7mm and 71 = 7c, whereas in the strong 
cooling regime, s = 2, 70 = 7c and 71 = 7mm- This method for deriving the 
behaviors of the breaks in the electron and synchrotron spectrum for the simple 
synchrotron blast-wave model were originally given by Sari, Piran, and Narayan 
(1998). [Exercise. Derive the temporal and spectral behaviors for the elementary 
blast-wave model consisting of an impulsive adiabatic blast wave with a strong 
forward shock sweeping a uniform surrounding medium.] 

The synchrotron shock blast-wave model has been used to fit afterglow 
data and deduce microphysical and environmental parameters. Detailed leptonic 
models include a synchrotron self-Compton (SSC) component, which is highly 
sensitive to the baryon- loading parameter Fq. 

3.4. Beaming Breaks and Jets 

An observer will receive most emission from those portions of a GRB blast 
wave that are within the Doppler angle Ob ~ 1/F to the direction to the ob- 
server. As the blast wave decelerates by sweeping up material from the external 
medium, a break in the light curve will occur when the jet opening half-angle 
6j < l/F. This is due to a change from a spherical blast wave geometry to 
a geometry defined by a localized emission region). Assuming that the blast 
wave decelerates adiabatically in a uniform surrounding medium, the condition 
ej ^ l/r = To\xbr/xdf/^ = T^\hr/tdf'^ implies 

tbr ~ 45(1 + z) (:^)V3^8/3 ^^yg ^ ^24) 
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from which the jet angle 



O.li^f^ i^f^ (25) 



can be derived. Note that the beaming angle is only weakly dependent on no 
and Eq. 

Numerical models show X-ray beaming breaks hidden by the effects of the 
SSC component. The important discovery (Frail et al. 2001) of a clustering, 
beaming-corrected energy for LGRBs opens the possibility to perform cosmo- 
logical studies with GRB data. 

3.5. Relativistic Shock Physics 

The structure of a shock is determined by continuity of the particle number, 
energy and momentum fluxes across the shock front. The pressure p = (7 — l)efce, 
where Cke is the kinetic energy density of the fluid, and 7 corresponds to the ratio 
of specific heats. For a nonrelativistic monatomic ideal gas, 7 = 5/3, whereas 7 
= 4/3 for a relativistic gas. For strong shocks, n' = 4rno. 

The equality of kinetic-energy densities at the contact discontinuity implies, 
for fluids made primarily of proton-electron plasma, that 

^ nfCT - 1) ^ nAT - 1) ^ AnoT^ ^ 4n(x)(T^ - F), (26) 



where the reverse shock Lorentz factor is F = y 1 — /^^ . The relativistic shock 
jump conditions for an isotropic explosion in a uniform surrounding medium 
imply (Sari & Piran 1995) 



= F = 




NRS 
RRS 



no Airx'^Tlnompc'^ A x^F^A 

""^(27) 

The relations between the forward shock (FS) and reverse shock (RS) Lorentz 
factors can be derived in the limit of a nonrelativistic reverse shock (NRS) and 
strong forward shock, and in the limit of a relativistic reverse shock (RRS)_and 
relativistic forward shock. The RS power is dE'/dt'\Rs = A{x)n{x)mpC' I3{V'^ - 
F). With the shock jump condition, one finds that {dE' /dt'\Rs)/{dE /dt'lps) = 
P, so that roughly equal power is dissipated as internal energy in the forward 
and reverse shock during the RRS phase. 

3.6. External Shock Model 

This is sufficient blast-wave physics that evolving forward and reverse shock 
emissions can be calculated for comparison with data. Synthetic light curves 
in the external-shock model can be derived and used to make predictions for 
the generic behavior of GRBs with smooth light curves. Spikiness of the light 
curve could originate from inhomogcncities in the external medium, before the 
blast wave has entered the adiabatic deceleration regime. Consider a blast wave 
intercepting a cloud with size r <g R/Tq that is located at an angle with 
respect to the line of sight to the observer. The duration of the received pulse 



20 



Dermer Fryer 



of radiation depends on the light travel-time delays from different portions of 
the blast wave as it interacts with the cloud. Photons emitted when the blast 
wave passes through the near and far sides of the cloud are received over a radial 
timescale 

2r r 

= - ffc ■ ^^^^ 

Photons emitted from points defining the greatest angular extent of the cloud 
are received over an angular timescale 

^rO 

tang — • (•^"^j 

Note that if r — s- R/Tq and 9 — > I/Fq, then tang —>■ R/^qC, as expected. When 

6 fa l/Po, tang ~ ^otr ^ Except for those few clouds with 6 ^ I/Pq lying 
almost exactly along the line of sight to the observer, tang ^ Highly variable 
light curves with reasonable ( ^ 10%) efficiency can be produced in an external 
shock model (Dermer & Mitman 1999). 

In the external shock model, a single relativistic wave of particles interacts 
with inhomogeneities in the surrounding medium to accelerate particles that 
radiate the prompt 7 rays. A central requirement for strong radiative efficiency 
in an external shock model for the prompt phase is that a strong forward shock 
is formed; otherwise the Lorentz factor P <C Pq and the radiation is strongly 
debeamed. A strong forward shock is formed when the comoving shell density 
n{x) 3> Po^^o- 

3.7. Internal Shock/Colliding Shell Model 

In the internal shock model, an active central engine eject waves of relativistic 
plasma that overtake and collide to form shocks. The shocks accelerate nonther- 
mal particles that radiate high-energy photons. The relative Lorentz factor of 
the two shells with Lorentz factors Pi and P2, with P2 = C^i, C ^ 1) is 

P,eZ = PlP2(l-/3l/32) \{C + C')- (30) 

For mildly relativistic internal shocks with a range of relative Lorentz factors 
1 ^ ^rei, C ~ 10) the Lorentz factor P of the shocked fluid with adiabatic index 

7 = 5/3 in the explosion frame is 

T,f ^ 2Pi ^relF^^^ ^ V^lTl = , 

where F = n2/ni is the ratio of proper frame densities of shell 2 to shell 1 
when they intercept each other (eq. [27]), and the final expression assumes that 
722 ~ ni and Trei ^ 1- The proper shocked fluid number and energy densities 
are 

Hsf = (4P^e/ + 3)ni w 4:Treini , and e^/ = TreiUsf rripC^ . 

Now consider the elastic collisions of shell 2 with mass m2 intercepting shell 
1 with mass mi. In an elastic collision, the Lorentz factor of the merged shell is 



r„ ^ . . (31) 

y m2/P2 + mi/Pi 
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The efficiency to convert the directed kinetic energy of the shells into internal 
energy is 

(mi + m2)Tm , . 

The efficiency is greatest when the shells have comparable mass and r2 ^ Ti; 
otherwise 77 ~ few %. When the contrast between the F factors of the shells is 
large, r/ 10 - 20% is possible. 

The rapid X-ray decline can be explained by high-latitude emission after 
the central engine has been turned off. Writing the flux density Fi, oc u'^t^ gives 
the curvature relation, which assumes that the spectral shape of the radiated 
flux is constant within the Doppler cone. The curvature relation is a = 2 — /?. 

In the standard model for LGRBs with relativistic winds and colliding shells, 
X-ray flares are made when the GRB engine is restarted. Long-lasting GRB 
central engines can also involve continual injection scenarios with pulsars. The 
generic long-duration GRB light curve at ~ 1 keV, from Swift observations 
(Zhang et al. 2006; Nousck et al. 2006), can be divided into a prompt phase 
(0), a decay phase (I), a plateau phase (II), an afterglow phase (HI), and a jet 
break phase (IV), in addition to X-ray flares (V). Kinematic shapes of GRB 
pulses can be calculated for illuminated shells when the thickness and duration 
of illumination of the shell are varied. 

3.8. Leptonic GRB Physics: Summary 

The success of the fireball/relativistic blast wave model arises from its ability to 
explain 

1. Large energy releases in short times; 

2. Escape of 7 rays; 

3. Afterglows at various wavelengths from radio through 7 rays. 

The physics is widely applicable to many nonthermal systems, including blazars, 
microquasars, pulsar winds,. . . 

This hardly exhausts leptonic blast-wave physics. Other interesting physics 
involves thermal photospheres, the "line of death" and synchrotron jitter radia- 
tion, and the origin of the Amati and Ghirlanda relations. 



4. Hadronic Processes and Cosmic Rays in GRBs 

Acceleration of ultra-relativistic protons and ions is favored in a blast wave 
physics scenario at least as much energetic ions, insofar as (1 — eg) of the power 
in the dissipation region emerges in the form of magnetic field energy or ions 
and, if ions, with energy ^ Fq GeV. The outline for this lecture is 

1. Ultra-high Energy Cosmic Rays (UHECRs) 

2. Photohadronic Processes, Energetics, and Power 

3. Hadronic Blast Wave Theory 
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4. Cosmic Rays from GRBs 

5. Neutrinos from GRBs 

6. GRBs in the Milky Way 

4.1. Ultra-High Energy Cosmic Rays 

The year 1912 is a landmark date in space science when the cosmic radiation, 
a "penetrating radiation from above," was discovered by Victor Hess by flying 
electroscopes on a balloon. The cosmic-ray energy density at GeV/nucleon en- 
ergies is ucR ~ 10^^^ ergs/cm^, with the total cosmic-ray kinetic energy density 
modulated by the outflowing Solar wind on 22-yr Solar cycle. The knee of the 
spectrum is at ~ 3 PeV, the second knee is at particle energy E ~ 10^^'^ cV, 
the ankle (or dip) at E ^ 10^^-^ eV, and the GZK cutoff is at Eqzk ~ lO^^'^ 
eV. The value of ucr ~ 10~^^ ergs/cm^ at ^ ^ Eqzk- 

The Pierre Auger Observatory (PAO), located in the Mendoza Province in 
Argentina at « 36° S latitude determines the arrival directions and energies of 
UHECRs using a hybrid technique consisting of four telescope arrays to measure 
Ni air fluorescence and 1600 surface detectors spaced 1.5 km apart to measure 
muons formed in cosmic-ray induced showers. Event reconstruction using the 
hybrid technique gives arrival directions better than 1°, and energy uncertainties 
better than ^ 20%. Two important discoveries were made in 2007, namely 

• GZK cutoff with the HiRes Observatory (Abbasi et al. 2008) and the PAO 
(2008). 

• Clustering of arrival directions toward AGN in the supergalactic plane, 
with the PAO (2007). 

• Interesting though ambiguous results on composition were also announced 
by the PAO in 2007 (Unger et al. 2007). 

The GZK cutoff, now seen clearly with HiRes and the PAO, contrary to ear- 
lier AGASA results, is a consequence of exponentially increasing energy losses 
of UHECR protons due to photopion-producing reactions with photons of the 
CMBR. For UHECR Fe, strong photo-dissociation though, primarily, one and 
two-nucleon losses in giant dipolc resonance reactions, also produce a GZK cut- 
off. There is no strong evidence for a lighter composition at E ^ 10^^ eV in- 
ferred from data taken with the PAO to energies E ^ 4 x 10^^ eV, contrary to 
pre- Auger suggestions (Watson 2006). Extrapolation of particle physics uncer- 
tainties to large values of total CM energy make deductions about composition 
from shower data uncertain. 

A major result announced by the PAO collaboration in 2007 was the cor- 
relation of arrival directions oi E ^ QO EeV UHECRs with nearby, d ^ 75 - 100 
Mpc, AGNs in the Veron-Cetty & Veron (2006) catalog. A marked excess in 
the direction of Cen A has generated much interest in the possibility that Cen A 
is an UHECR source. Both radio-loud AGN and GRBs remain plausible candi- 
dates, but the absence of strong radio galaxies within !=a 100 Mpc of a number of 
UHECR arrival directions is puzzling for a radio-galaxy origin of the UHECRs. 



La Plata Lectures on GRBs and Fermi/GLAST 



23 



4.2. Photohadronic Processes, Energetics, and Power 

Like filling a bathtub, filling the Galaxy with cosmic rays, or the universe with 
UHECRs, is a balance between injection and escape (loss). Symbolically, 

ucr{E) ~ tioss{E)ecR{E) . (33) 

The emissivity ecR{E) is source dependent, and expected to favor strong non- 
thermal sources. Relevant photohadronic proceesses are photomeson and pho- 
topair production, and photodisintegration for ions. The photomeson or pho- 
topion cross section resembles a threshold step-function due to the onset of 
various resonant and nonresonant channels above threshold. Photopion losses 
due to interactions with photons of the extragalactic background light (EBL), 
importantly consisting of the CMBR, and Bethe-Heitler (A^7 Ne^e~) pho- 
topair losses on nucleon A^. In addition, one must consider ion-synchrotron and 
universal expansion losses. (Exercise: Derive the ion-synchrotron energy loss 
rate in a tangled magnetic field, and the maximum synchrotron photon energy 
and particle energy in Fermi acceleration scenarios.) For photopion losses, we 
use a step-function approximation, from which the energy loss mean-free-path 
for particles of energy E can be derived. 

The UHECR emissivity (or luminosity density), from eq. ([33]) . is ecR{E) ~ 
ucr{E) /tioss{E) . Using PAO measurements and results of energy- loss mean- 
free-paths, then ecR{E) ~ 10^^ ergs/Mpc^-yr for E ^ 10^^ eV. This is within 
an order of magnitude of the electromagnetic emissivity of GRBs, as noted in 
1995 by Vietri and Waxman. If the baryon-loading factor /b ^ 1, as implied 
when Ce <C 1, then LGRBs could power the UHECRs; SGRBs do not seem to 
have the required emissivity. 

LLGRBs, which take place at low redshifts (GRB 980425/SN 1998bw was 
at d ~ 40 Mpc), can also, in principle have the requisite emissivity. Estimates 
show emissivities of the outflowing kinetic energy in LLGRBs of ~ 250 x 10^^ 
ergs/Mpc^-yr (Wang et al. 2007); however only a small fraction of this luminosity 
density can be expected to emerge in the form of hard X-rays and soft 7 rays. 

4.3. Hadronic Blast Wave Theory 

The elements of blast- wave theory applied to leptons is, with appropriate changes, 
directly applied to hadrons. In the simplest model, cosmic-ray protons are in- 
jected downstream of the shock with spectrum 

7V'(7')cx7'-^^(7';7™n,7™aJ, (34) 

normalized to the number of swept-up protons A^'o = 47rx'^/3, the swept-up 
power, and the magnetic field B defined in terms of shocked fluid energy density 
(Lecture 3). The minimum comoving proton Lorentz factor 7rnin ~ ^(^^ 

), and 

the maximum is set by equating size scale with Larmor radius. For equipartition 
magnetic fields, cosmic-ray protons reach ultra-high, ^ EeV energies after freely 
escaping from the blast wave into interstellar space. [Exercise: work out these 
relations.) 

The emission spectrum from a GRB includes in addition to a leptonic com- 
ponent a photohadronic component. The three dominant collision processes are 
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(i) secondary nuclear production; (ii) photomeson production; and (iii) pho- 
topair production. Energetics arguments favor photohadronic processes over 
the nuclear collision processes. Proton and ion synchrotron radiation must 
also be considered. The cascade 7-ray spectrum is initiated by decay of sec- 
ondary mesons, attenuation of high-energy photons, and synchrotron radiation. 
Hadronic emission decays more slowly than leptonic emissions in standard blast- 
wave model calculations (Bottcher & Dermer 1998). {Problem: Analytically 
examine correlations in variability behavior for proton synchrotron 7 rays and 
leptonic synchrotron emission. Do the same for photohadronic processes.) 

4.4. UHECRs from GRBs 

The astrophysics to calculate the UHECR energy spectrum measured here at 
Earth can be extensive (Bcrczinskii & Grigor'cva 1988). Assuming a homoge- 
neous source type that explodes and releases the same UHECR spectrum into 
intergalactic space throughout cosmic time enormously simplifies the problem. 
The star formation rate factor for GRB progenitors could follow the classical SFR 
(Hopkins & Beacom 2006), or have a different dependence due, e.g., metallicity 
effects. For very active SFRs at z ^ 1 compared to now, the pair-production 
trough becomes more pronounced to explain the ankle feature through photopair 
losses. The model of Wick, Dermer, and Atoyan (2003) for the UHECR spec- 
trum from GRB sources predicts a very sharp GZK cutoff but requires a large 
baryon load, « 30 - 100. (Note that the lower normalization of Auger vs. 
HiRes brings the value of down.) 

Different choices for the GRB star formation rate, which normalizes the 
luminosity density (emissivity) at different z, can be normalized to various SFR 
factors. Very active star formation rates can be ruled out by comparing Swift 
and pre-Swift statistical GRB data (Le & Dermer 2007). The UHECR spectrum 
is also quite sensitive to the maximum energy of accelerated ions or protons. 

4.5. Neutrinos from GRBs 

The opening of the high-energy (> TeV) neutrino window is anxiously awaited, 
when statistically significant detection of a cosmic point or extended source 
of high-energy neutrinos with a deep ice-based or water-based Cherenkov light 
sensitive detectors happens. Charged-current interactions of a neutrino with a 
nucleon induce muon, electron, or tau production(Gaisser et al. 1995). These 
particles cascade and shower to make Cherenkov light. Detector optical modules 
(DOMs) in the IceCube detector at the South Pole detect upward-going tracks 
to screen out the intense background from downward-going showers induced by 
cosmic-ray muons. 

Significant extraction of UHECR energy via photohadronic processes can 
be made in a collapsar- model GRB, with a significant fraction going into neutri- 
nos. The efficiency for neutrino production depends importantly on the baryon 
loading and Doppler factor. A spectrum with index —2 minimizes energy require- 
ments and gives detectable neutrino production for the large GRB baryon load- 
ing required to fit the UHECR spectrum. Interesting anti-correlations between 
neutrino and 7-ray brightness arise because of the underlying physics. Photome- 
son cascade calculations consist of many generations of Compton-upscattered 
and lepton synchrotron radiation. A 7-ray spectral model of a GRB must in- 
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elude, at least, a photohadronically induced 7-ray component from UHECRs in 
GRB blast waves, in addition to the lepton synchrotron and SSC radiation. The 
high-energy neutrinos are formed between ~ 100 TeV and 100 PeV. 

Photomcson production in a GRB blast wave, for parameters used to fit the 
UHECR spectrum, leads inevitably to a neutral beam of neutrons, 7 rays, and 
neutrinos. Many interesting implications of the neutral beam model (Atoyan & 
Dermer 2003) follow, first being detectability of GRBs in high-energy neutrinos 
for reasonable collapsar model GRB parameters, though with Doppler factors 
;S 100. Subsequent photopion interactions from UHECR neutrons or neutron- 
decay UHECR protons induce beam of 7 rays and leptons that cascade, making 
hyper-relativistic, highly polarized synchrotron radiation. A classical LGRB 
model for UHECRs allows fairly definite predictions to be made involving only 
a few parameters for the cosmogenic GZK neutrino spectrum. This gives the 
"guaranteed" UHECR u spectrum for a given astrophysical/cosmological model 
of UHECR origin. 

4.6. GRBs in the Galaxy 

Cosmic Rays from GRBs in the Galaxy. Depending on the beaming factor 
(pbmy the rate of GRBs in Milky Way is estimated to be ~ 1 per (0.1 - 1) Myrs. 
The rate of GRBs producing astrobiological effects at Earth is ~ 1/Gyr. 

For a beaming factor ~ 1/50 1/500 (Frail et al. 2001) the mean 7-ray 
energy in X/7-rays is ~ 5 x lO^'' ergs. The likelihood of a recent ( ^ Myr) 
GRB in our Galaxy is scaled from BATSE rate of f» 2 GRB/day. From these 
estimates, !=a 0.3 - 1% of SNe collapse into black holes, implying ~ 1 GRB every 
~ 3 - 100 kyrs in the Galaxy The expected number of recent GRBs within 
r(kpc) of Earth over a period of duration t is, roughly, 

{Ngrb) ^r\kpc)t{Myi) . (35) 

{Exercise: Verify, correct, or improve these estimates.) 

Upon injection by a GRB, high-energy cosmic rays diffuse throughout the 
disk into the halo of the Milky Way. For isotropic turbulence and a diffusion 
coefficient D{E), a particle cloud diffuses according to the relation 

N{E, r, t) oc r^,J^ exp (-r^/rl/y) , r^iffiE) = ^j2D{E)t. (36) 

A two-component turbulence spectrum, corresponding to turbulent energy injec- 
tion at the pc and 100 pc scales following Kolmogorov and Kraichnan behaviors, 
was used to model (Wick et al. 2004) propagation of cosmic rays and ions be- 
tween K, 100 TcV and ?a lO^'' cV. The diffusion coefficient is limited by the 
relation rdiff{E)/t < c and, in general, is anisotropic. 

In its simplest form, such a model is sufficiently quantitative to fit the 
energy spectra of UHECR ions, as measured with the Karlsruhe observatory, 
KASCADE (Antoni et al. 2004). The injection spectrum of the ion component 
is fixed, with composition varied to improve the fit. The propagation character- 
istics depends of the rigidity coefficient, essentially energy per charge at these 
energies. A crucial issue is anisotropy. The combined Galactic GRB model 
for cosmic rays through the knee, combined with the extragalactic/cosmological 
UHECRs make a complete model for high-energy cosmic rays. Propagation of 
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high-energy cosmic rays in the Galaxy exhibit surprising but easily understood 
effects. 

Astrobiological Effects from GRBs in the Galaxy. The rate of intense events 
from GRBs in the Galaxy can be estimated from the previous results (Dermer 
&: Holmes 2005). Define the bolometric photon fluence (p = SipQ with reference 
to the Solar energy fluence ifQ = 1.4 x 10^5 ergs cm~^ received at Earth in 
one second. Significant effects on atmospheric chemistry through formation of 
nitrous oxide compounds and depletion of the ozone layer is found when S ^ 10^ 
- 10^ because of the very hard incident radiation spectrum of GRBs that is 
reprocessed into biologically effective 200 - 320 nm UV radiation. 

Using the standard energy reservoir for LGRBs to establish apparent energy 
release for a jet with opening half-angle 9j, one finds that the maximum sampling 
distance Rs of a GRB with apparent isotropic 7-ray energy release E^^iso to be 
detected at the fluence level ip > iptu = is 



„ _ l E^,iso ^ 1-1 kpc / £51 

^ "i^n^th' i9,/0.1)^{S/10^)- ^^'^ 

If one GRB occurs every 10^t5 years in the Milky Way, the rate of biologi- 
cally significant events is 

Thus a GRB at a distance « 1 kpc with <S S> 10^ takes place about once every 
Gyr, and more frequently if ts = 0.1. 

A GRB pointed towards Earth produced a lethal flux of high-energy photons 
and muons that destroyed the ozone layer, killed plankton, and led to trilobite 
extinction in the Ordovician Epoch (Melott et al. 2004). Geological evidence 
points toward two pulses: a prompt extinction and an extended ice age. A 
muon dose and intense flux of ionizing photons from a GRB could have produced 
the prompt extinction. Delayed cosmic rays could have produced the later ice 
age. {Exercise: Calculate muon dose at Earth's surface from UHECR neutron 
impacts on the upper atmosphere made by a GRB in the Galaxy. Describe 
effects.) 

The issue of metallicity-dependence of LGRB progenitors introduces a large 
additional uncertainty into these estimates (Stanek et al. 2006). Local LLGRBs 
are found in low-metallicity host galaxies. This may also be the case for LGRBs. 
Recent work has also considered the role of LLGRBs as sources of UHECRs. 

A complete model for cosmic rays developed by Wick, Dermer, and Atoyan 
(2004) has 

1. Cosmic Rays below ^ 10^^ eV from SNe that collapse to neutron stars 

2. Cosmic Rays above ^ 10^^ eV from SNe that collapse to black holes 



3. CRs between knee and ankle/second knee from GRBs in Galaxy 
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The highest energy cosmic rays originate from outside our Galaxy, because their 
Larmor radii exceed the size scale of the Galaxy. (Problem: Derive the transition 

energy from galactic to cosmological dominance of cosmic rays for a realistic 
galactic magnetic field model. Fit to data using simplifying assumptions.) 

4.7. Hadronic GRBs: Summary 

The acceleration of ultrarelativistic hadrons in GRB blast waves introduces 
many new aspects to the GRB problem. Specializing to LGRBs and UHECR 
protons, we find that 

• LGRBs are a viable source of UHECRs; 

• Hard 7-ray emission components are formed by hadronic cascade radiation 
inside GRB blast wave; 

• A second emission component is formed by outflowing high-energy neutral 

beam of neutrons, 7-rays and neutrinos; 

• GRBs can be detectable high-energy neutrino sources, which would con- 
firm UHECR acceleration; 

• Cosmic rays can be accelerated by GRBs in the Galaxy; and 

• GRBs could be responsible for ionizing or extinction events on Earth. 

A Northern hemisphere Auger, like its Argentine counterpart, will produce the 
first full-sky map of the universe not in the electromagnetic window, but in the 
particle window. 

5. GRB and 7-ray studies with Fermi/GLAST 

Each major advance in GRB science is a result of new instrumentation. The 
next GRB epoch has already began with the launch of GLAST on June 11, 2008, 
renamed the Fermi Gamma-ray Space Telescope (FGST) on August 26, 2008. 
Besides GRBs, the FGST will vastly increase our knowledge of astronomical 
7-ray sources. The historical perspective presented here will soon be superseded 
by the FGST, but gives us the opportunity to guess what it might see. 

Besides the steady diffuse glow of the Milky Way at ICQ MeV - GeV energies 
from the decay of pions formed as secondaries in cosmic-ray collisions with dust 
and gas, the 7-ray sky is pulsings from pulsars, flarings from the Sun and blazars, 
burstings from GRBs, and revealing 7-ray enhancements from sources in the 
Solar system. Galaxy, and beyond the Galaxy at sites of cosmic-ray production 
or dark-matter annihilation. 

The FGST will join an ongoing revolution in high-energy astronomy, in- 
cluding ground-based air and water Chercnkov telescopes VERITAS, HESS, 
Cangaroo HI, MAGIC-2, and Milagro and its successor HAWC. At GeV ener- 
gies, AGILE (with its super- AGILE X-ray detector), will be joined by the FGST 
with its GBM and LAT. This is in addition to multimessenger observatories (in- 
cluding Auger, KMSNeT, HiRes, LIGO, LISA, ANITA) and X-ray telescopes 
(e.g., Chandra, XMM, RXTE, Suzaku ) {Exercise: Report on various obser- 
vatories.) 
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5.1. EGRET and Fermi/GLAST 

The Energetic Gamma Ray Experiment Telescope (EGRET) on CGRO was 
most sensitive between ~ 100 MeV - 5 GeV with a field-of-view equal to about 
1/24*^^ of the full sky. Its maximum effective area between 100 MeV and 1 
GeV was about 1200 cm^, and its point-spread function was 5.7° at 100 
MeV, improving oc E^^^"^ at higher energies. For a nominal two-week observing 
period, EGRET reached integral fluxes ^ 15 x 10"® ph(> 100 MeV) cm"^ s~^ 
for liigh-latitude point sources. The FOV is defined as the subtended solid angle 
swept from the zenith to the angle where the effective area is one-half of on-axis 
effective 

EGRET provided tantalizing clues to 7-ray emission from GRBs, including 
GRB 940217 with a tail of 100 McV - GcV emission to ~ 93 minutes after the 
GRB trigger (Hurley ct al. 1994). A 20 GcV photon appeared after the end 
of Earth occultation. Forming crude spectra during prompt phase and after- 
glow phase shows clear evidence for a second, high energy emission component. 
Moreover, ~ 100 MeV 7 rays are found during peak pulses. 

Not only GRB 940217, but the superbowl GRB 930131, which could be a 
SGRB, displayed hard emission after the end of the prompt phase. The proto- 
type Milagro experiment presented evidence for multi-TcV emission from GRB 
9704171A, which can best be confirmed if new members of this class were de- 
tected. 

Most important for expectations with the FGST is to consider results for five 
GRBs from the EGRET Spark Chamber. The average spectrum of high energy 
GRB emission within 200 seconds of BATSE trigger for GRB 940301, GRB 
940217, GRB 930131, GRB 910601, GRB 940503 is well fit with a differential 
power law photon spectrum with index 1.951q;4^ (Dingus et al. 1998). Seven 
GRBs were detected with EGRET either during the prompt sub-MeV burst, 
or after sub-McV emission has decayed away. The ratios of the 100 MeV - 5 
GeV EGRET fluence to > 20 keV BATSE fluence are typically a few percent 
to a few tens of percent. (Note that these GRBs were all observed early in 
the EGRET mission before the spark chamber gas was degraded. Silicon strip 
detector technology in the FGST avoids the use of consumables.) 

The Total Absorption Calorimeter on EGRET scintillated in response to 
very bright GRB events, even for off-axis GRBs. Joint BATSE/EGRET TASC 
analysis discovered anomalous 7-ray emission components in GRB 941017 not 
easily explained with the leptonic blast- wave model (Gonzalez et al. 2003). 

The FGST, as already noted, consists of two telescopes, the Large Area 
Telescope, the LAT, and the Fermi GLAST Burst Monitor, GBM. The sensi- 
tive energy range of the LAT is between 50 MeV and 100 GeV (self-vetoing 
backsplash limited the highest energies of EGRET to 5 GeV), with a FOV of 
1/5*'^ of the full sky, and a PSF of ^ 3.5° (0.55°) at 100 MeV (1 GeV). Its 
on-axis effective area for GeV photons is « 9000 cm~^, and the nominal ob- 
serving strategy is to scan the sky every three hours. The combined effective 
area and smaller FOV means that the FGST can reproduce EGRET's one-year 
capabilities within about 4 days, and will reach one-year detection thresholds of 
^ 0.4 X 10~® ph(> 100 MeV)/cm2-s. The GBM can, to first order, be consid- 
ered comparable to a slightly smaller BATSE, detecting 200 GRBs per year. 
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but with better sensitivity between « 1 - 30 MeV due to the BGO (Bismuth 
Germanate) scintillator, as well as sensitive down to ~ 8 kcV. 

The effective area for the the FGST telescope depends on the observing 
mode, bTit naive comparison of effective areas and point spread functions reveals 
how superior the FGST will be with respect to EGRET. (For the latest Fermi 
LAT instrument performance, google "GLAST LAT Performance.") Thin and 
thick sections refer to different thicknesses of conversion layers in the LAT tracker 
to optimize for effective area (thick layers) or direction (thin layers). Different 
predictions about FGST LAT detection of GRBs can be made by scaling from 
the relative ratios of fluences in BATSE to EGRET for the 5 spark chamber 
GRBs. About 20 GRBs full sky per year with more than 10 (> 100 MeV) 7 
rays, or about four per year in the LAT FOV with more than 5 (> 100 McV) 
7 rays (excluding autonomous slewing maneuvers) , are predicted (Le &; Dermer 
2008). Only about 1 GRB with more than 100 (> 100 MeV) 7 rays is predicted 
through the five-year nominal Fermi/GLAST lifetime, but presence of second 
components or new classes of GRBs should improve detection rate. 

FGST data is proprietary until one year after the start of Phase I, the 
first year of science operations. Even during the first year, the GLAST team is 
obligated to release light curves and spectral data as soon as practical for 7-ray 
transients that exceed a flux of 200 x 10~^ ph(> 100 MeV)/cm^-s, which 
should occur at the rate of once every week or so. On top of that, the light 
curves and spectral behavior of some 24 sources of interest will be released to 
the community. 

Analysis of EGRET data to search for point and extended points of radi- 
ation requires, especially at low (|6| < 10°) Galactic latitudes, a diffuse model 
for cosmic-ray interactions in the Galaxy. The third EGRET catalog (Hartman 
et al. 1999) lists 271 7-ray sources, including the single 1991 solar flare bright 
enough to be detected as a 7-ray source, the Large Magellanic Cloud, five pulsars, 
one probable radio galaxy detection (Ccn A), and 66 high-confidence identifi- 
cations of blazars (BL Lac objects, fiat-spectrum radio quasars, or unidentified 
flat-spectrum radio sources). In addition, 27 lower confidence potential blazar 
identifications are noted. Finally, the catalog contains 170 sources not yet iden- 
tified firmly with known objects. 

Surprises and important results are expected from FGST regarding 

1. 7 rays from associations of high mass stars with strong stellar winds; 

2. Spectra of 7-ray emission from supernova remnants, allowing 7-ray astron- 
omy a flnal chance to answer the question of Galactic cosmic-ray origin; 

3. Origin of microquasar 7-ray emissions, whether due to a scaled-down mi- 
croquasar jet or to particles accelerated the shocks formed at the interface 
of a binary system consisting of a young pulsar and a high-mass star; 

4. Galactic astronomy, using other data bases (e.g., WMAP) to search for 
gaseous components in the Galaxy that are illuminated by cosmic-ray in- 
teractions; 



5. 



Search for 7 rays from normal, starburst, and IR luminous galaxies; 
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6. 7 rays from cosmic rays energized by structure formation shocks in clusters 
of galaxies; 

7. 7 rays from dark matter annihilation; 

8. Blazars and their misaligned radio galaxy counterparts. 

5.2. AGN Studies with the Fermi Gamma-ray Space Telescope 

The standard blazar model features collimated ejection of relativistic plasma 
from supermassive black holes. The relativistic motion accounts for lack of 
77 attenuation, as in the case of GRBs. Other evidence for relativistic outflows 
include superluminal motion, super-Eddington luminosities, high-energy beamed 
7 rays made in Compton or photo-hadronic processes. There is a considerable 
difference in the environments of GRBs and blazars, for example, the intense 
external radiation field from broad line-region gas in FSRQs. 

The redshift distribution of EGRET blazars leads to model predictions for 
the FGST. Connections between different types of supermassive black holes with 
jets of radio and 7-ray emitting plasma can be made based on data with the 
FGST, including tests of the blazar main sequence, and studies of black-hole 
engine and jet physics. FGST observations with correlated multiwavelength 
campaigns will provide important data for modeling studies, and for searching 
for anomalous blazar 7-ray emission components, including orphan flares and 
hadronic 7-ray signatures. Long straight radio and X-ray jets, like in the FR2 
radio galaxy Pictor A, could reveal UHECR acceleration in blazars through 
neutral beam processes. Acceleration of UHECRs in blazars, like in GRBs, will 
be most decisively demonstrated with neutrino detection. The detection of one 
or two ~PeV neutrinos from a blazar during flaring conditions will overturn our 
thinking about how radio lobes are formed, because a decaying neutrons and 
attenuated UHE 7-rays can power the knots, hot spots, and lobes of a radio 
galaxy. Detection of pair halos from nearby radio galaxies would give evidence 
in favor of the UHECRs made in GRBs. 

The superposition of 7-rays formed by various sources throughout cosmic 
time produce an unresolved 7-ray background intensity. Only about 10 - 20% 
of the diffuse background can be made by FSRQ and BL Lac blazars, based on 
analysis of the EGRET results (Dermer 2007). Other source classes that can 
make up the diffuse extragalactic 7-ray background include 

1. Star-forming galaxies; 

2. Starburst galaxies; 

3. Pulsars; 

4. Galaxy cluster shocks; and 

5. Dark matter annihilation. 

Although GRBs are momentarily very bright, they are so rare that their con- 
tribution to the diffuse background is small. (Problem: Construct fully analytic 
or semi-analytic [solution in quadrature, i.e., single integral] models of various 
source classes.) 
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5.3. Concluding Remarks 

Some important problems in GeV 7-ray astronomy that will soon be opened for 
study in view of anticipated FGST discoveries include 

1. Particle acceleration theory; 

2. Origin of the galactic cosmic rays; 

3. Jet physics, and the differences between radio/7-ray black-hole sources; 

4. Blazar demographics; 

5. Search for hadronic 7-ray emission components and the sources of UHE- 
CRs; and the 

6. Origin of the diffuse/unresolved 7-ray background. 

In view of the Fermi Gamma ray Space Telescope for GRB studies, we can 
expect a new golden era. 
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